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We study the absorption of a massless scalar field by a static black hole. Using the continuity 
equation that arises from the Klein-Gordon equation, it is possible to define a normalized absorption 
rate r(f) for the scalar field as it falls into the black hole. It is found that the absorption mainly 
depends upon the characteristics wavelengths involved in the physical system: the mean wavenumber 
and the width of the wave packet, but that it is insensitive to the scalar field's strength. By taking 
a limiting procedure, we determine the minimum absorption fraction of the scalar field's mass by 
the black hole, which is around 50%. 



o 

1—5 



o 

I 



> 

m 

cn 

o 



X 



PACS numbers: 04.40.-b,04.25.D-,95.35.+d,95.36. 



I. INTRODUCTION 

Black holes, a concept that emerged from the simplest 
exact solution of Einstein's equations, are some of the 
most fascinating objects in gravitational physics. Equally 
fascinating is our current belief that most galaxies must 
host a supermassive black hole (SMBH) in their cen- 
ter, with mass values in the range of 10^ to 10^" solar 
masses, most likely in a state of very low matter accre- 
tion nowadays [1,0]- In particular, the measurements of 
the velocities of stars near the center of the Milky Way 
have provided strong evidence for the presence of a SBH 
with a mass of around 4 x IO^MqQ- 

There are some models that attempt to explain the 
present existence of galactic SMBH's. Among others, we 
can mention the collision of two or more black holes to 
form a larger one, the core-collapse of a stellar cluster, 
and the formation of primordial black holes directly out 
from the primordial plasma in the first instants of time 
after the Big Bang[3i]. 

The key point in the discussion are the features of the 
precise mechanism under which a black hole can accrete 
enough matter to become supermassive. In particular, 
some authors have proposed that primordial black holes 
(PBH) can go supermassive simply by accreting matter 
from a cosmological scalar field related to dark energy 
(quintessence). In a first study, the authors in[4j (see 
alsop-Q) found that PBH could have effectively accreted 
enough matter from a quintessence field endowed with an 
exponential potential. 

The calculations for the accretion were in fact based 
upon the simple and exact results of the accretion of a 
massless scalar field into a black hole found inQ], see 
also [3, [a. [Toll . However, the results in[J| were later re- 
fused in[ll|, where was shown that the quintessence fiux 
must decrease slower than t~^ for PBHs to grow at all. 
This same result seems to have been confirmed by other 
authors under more general assumptions |9l.llCl Il2l). 
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On the other hand, a related topic is the use of a cos- 
mological scalar field as model for dark matter in the 
Universe [13[, and the possibility that they can be the 
dominant matter in galaxy halos|14|. If so, then one has 
to address the accretion of this dark matter scalar field 
into the central SBH that seems to be present in most 
galaxies |15llla |. 

The aim of this paper is to present some simple results 
of the interaction of a scalar field with a black hole, with 
numerical calculations based upon previous works in the 
literature [l7| that may be useful in the understanding 
of the accretion, in general terms, of cosmological scalar 
fields into black holes. 

We shall make use of the fact that there exists a con- 
tinuity equation of the scalar field as long as the back- 
ground spacetime is static 17[ • This fact will allow us to 
quantify the absorption rate of a scalar wave packet by a 
black hole in a more precise manner in terms of absorp- 
tion fiux and decay rates. For simplicity, we will only 
focus our attention in the case of a massless scalar field. 

A brief summary of the paper is as follows. In Sec. [Tl] 
we set the mathematical background for the equations of 
motion, boundary conditions, and initial conditions for 
the scalar field's wave packet. Here we also show the 
existence of a continuity equation arising directly from 
the equation of motion of the scalar field. In Sec. IHI[ we 
present the main numerical results, and the description 
of the fall of the scalar field in terms of a normalized ab- 
sorption rate. The latter arises naturally from the use of 
the continuity equation found in Sec. [Hi Finally, Sec. IIVI 
is devoted to conclusions and final comments. 



II. MATHEMATICAL BACKGROUND 

We first consider a fixed Schwarzschild background 
with an Eddintong-Finkelstein (EF) gauge, which is de- 
fined such that t -I- r is an ingoing null coordinate. Using 
the 3-t-l decomposition of the metric [17l[l8j . the 3-metric 

lij is 



Itj 



diag [a^(r), r'^,r'^ sin^ 9] 



(1) 



where a^(r) = 1 + 2GM/r, G is Newton's constant, and 
M denotes the mass of the black hole. The lapse a and 
shift /3* = [/3, 0, 0] functions are, respectively. 



a(r) = a ^{r) , (i{r) 



2GM 



(2) 



It is illustrative to calculate the coordinate velocities of 
null geodesies, that are given by 



dr 
'di 



c± 



~P±a/a. 



(3) 



Notice that the use of the EF gauge, from Eqs. ([2]), is 
manifest through the condition c_ = — 1 for all points in 
the background spacetime (we use units in which c—\). 
The Klein-Gordon (KG) equation for a massless self- 
interacting scalar field is 



-gg^^'^dM^Q. 



(4) 



In order to solve it, it proves convenient to define two 
first order variables |14l Il7| . 



n(t, r) = drcj) , 



^{t,r)^-{dt4>-m , (5) 

a 



with the help of which Eq. ^ is represented by the fol- 
lowing three first order equations 

- (9(0 - I3dr4>) = ?7 , (6a) 

a 

- {dtn - /39,n) = drfi + 1-7? + a,/3 n , (6b) 

a Zr 

-,2 



1 



(9(77 - (3drii) 
a a 



2~ 



1 n 



+ Kt] , (6c) 



where K is the trace of the extrinsic curvature. Eq. (|6a|) 
arises from the very definition of 77 and 11, whereas the 
equation for 11, Eq. (|6b[) . arises from the combination of 
Eqs. ([U; that of 77 in Eq. (|6cl) arises from the original KG 
equation ^. 

We shall take the quantity [GM) as the unit for dis- 
tance and time, so that the radial and time coordinates 
are made dimensionless through the change r — >■ {GM)f 
and t — >■ {GM)i, where a hat denotes dimensionless vari- 
ables. Notice that the unit for distance and time is half 
the usual Schwarzschild radius, rs = 2GM. The scalar 
field is made dimensionless by the change — >■ 4>/y/G. 
Accordingly, the rest of the scalar field variables should 
be changed by the expressions 11 — >■ {m?pi/M)Ii, and 
rj — ?► {mpi/M)?), where the Planck mass is defined as 



'Pi 



G- 



We will impose outgoing-radiation boundary condition 
upon our field variables at the outermost points of the 
numerical grid. As for the innermost points, as long as 
they are inside the event horizon, there is no need to put 
a boundary condition, because the light cones there point 
inwards in the EF gauge, that is, c+ < and c_ = —1, 
see Eqs. ©, and ©. 



The initial data for the scalar field in our numerical 
experiments will be a Gaussian profile modulated by a 
spherical wave of the form 



(^) ^ ^£^iiMe-(^-^o)V*^ 



(7) 



which is centered at f = fo, and has amplitude A and 
width a. The Gaussian distribution of wavenumbers in 
Fourier space has a mean value (k) = kg, and variance 
((/c — fco)^) = f/(T^. It should be noticed here that the 
wave number was made dimensionless by the change k — >■ 
k/{GM). 

On the other hand, by means of a lengthly but oth- 
erwise straightforward calculation, it can be shown that 
the KG equation Q can be written in the form of a con- 
tinuity equati' 



1 

ition (HI) 

tionflTl 



I 



dtp-^df{f'J' 



0, 



(8) 



where the charge density p and the scalar field current 
density J^ are, respectively. 



r 




(9b) 



In what follows, we will skip the hats of the variables, in 
the understanding that they have been made dimension- 
less. We may use a hat again on the variables whenever 
confusion may arise. 



III. NUMERICAL RESULTS 

The massless scalar field corresponds exactly to the so- 
lution of the (homegeneous) wave equation in a curved 
spacetime. One important feature of the massless case, 
which is helpful for the study of wave packets, is that the 
field keeps its shape as it falls into the black hole. Heuris- 
tically, this can be seen from the fact that the phase Vp 
and group Vg velocities of the wave packet are both equal 
to that of light, Vp — Vg = c±. For illustration purposes, 
this nice feature can be seen directly in the motion of the 
wave packets shown in Fig. [TJ 

We are solving differential equations with a finite dif- 
ferencing method, which involves the truncation of a Tay- 
lor series expansion. It is then necessary to show the 
proper convergence of numerical output as the spatial 
grid is refined. In Fig. [21 we show numerical runs with 
three different resolutions, where Ri has coarse resolu- 
tion, i?2 has medium resolution, and R3 is the finest. 
As expected, the runs show that the numerical code is 
second order convergent. 

To study the rate at which the a scalar field wave 
packet is absorbed by the black hole, we rely on the con- 
tinuity equation. Notice that Eq. ([S]) looks pretty much 
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FIG. 1. Motion of a Gaussian packet for different values of its 
mean wavenumber ko, see Eq. ((Tj). The wave packets retain 
their shape as they approach the black hole's horizon, located 
at r = 2. The wave packets were given ingoing initial condi- 
tions, 77 = n. The time in each row proceeds from right to 
left. 
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FIG. 2. Convergence plot for the numerical code used to solve 
the KG equation in a fixed Schwarzschild background. The 
example corresponds to the maximum of the wave function 
<l>max at each time in the simplest case fco = (see also Fig.[l}. 
Shown are three different runs with resolutions _Ri: Ar = 
0.025, R2: Ar = 0.0125, and Ry. Ar = 0.00625. The plots, 
after the indicated scaling, agree in their profile, then the 
numerical output is second order convergent. 



the same as a typical conservation equation in flat space- 
time. Taking a bounded proper volume V, we find 

dt f V^pdrdn^dt [ Anr^ pdr = 4:TT {r^.r)p . 
Jv Jv ^^ 

(10) 

Under the assumption that the scalar field current decays 

rapidly enough as r2 — >■ 00, we find the useful result 

Ma dt Ma ^ ^I'-S ^ ^' ^ ' 



where M^ = 47r Jy r"^ p dr is the total scalar field mass 
contained in the proper volume V. In fact, M^ is the 
conserved charge of the field as suggested by the conti- 
nuity equation ([5]). We can monitor the absorption rate 
of the (total) mass of the wave packet by calculating the 
scalar field current going through the inner surface; in our 
case, the inner radius is the black hole's horizon, ri = rs- 

Following standard notation in Physics, we have de- 
noted the decay rate of the wave packet's mass as r(i), 
whose units are given in terms of {GM)~^. In our case, 
this decay rate is just the normalized fiux at the horizon 
of the black hole, being the normalization factor the to- 
tal mass of the wave packet that still remains outside the 
black hole's horizon. 

Typical curves for the decay rate are shown in Fig. |31 
An interesting and unexpected result is that the (normal- 
ized) decay rate T{t) does not depend on the Gaussian's 
amplitude A, that is, it does not depend on the field's 
strength. This means that larger packets are absorbed 
at the same rate as are smaller packets. We can notice 
though that the mean wavenumber has an effect on the 
absorption, as the latter increases for larger values of kg. 

If we integrate Eq. ((TT|) . we can find the total mass 
outside the horizon as a function of time, 



M^t) 



M 



0,i 



exp 



T{t)dt] =exp(-A(i)) 



(12) 



where M^^i is the initial total mass, and A(i) is the (ex- 
ponential) absorption ratio. If the integral is calculated 
for the total time the wave packet is interacting with the 
black hole, it should give us the total absorption ratio of 
the wave packet. For the cases shown in Fig. [31 we have 
found that absorption is about 92% for fco = 0, and 100% 
for fco = 2. As a matter of fact, our numerical experi- 
ments showed that total absorption is always achieved if 
ko > rg^ (see alsojjj). 

As a final step, we study the dependence of A on the 
width of the wave packet; for definiteness, we focus our 
attention in the case fco = 0, which is also the most 
dispersive one. Numerical results are shown in Fig. 31 
and we notice that the absorption decreases as the wave 
packet becomes wider. It can be verified that the points 
can be fitted by a function of the form 



A(cr) = Cqc' 



62: 



(13) 



where, in the present case, a fitting procedure shows that 
eo = 4.98, ei = 0.259, and 62 = 0.67. In particular, if 
wave packets as wide as necessary were allowed, then 
Eq. p^ suggests that the total absorption would be 



lim A(cr) = 0.67. 



IV. CONCLUSIONS 



(14) 



The motion of a wave packet in a black hole space- 
time raised some interest in the cosmological community 
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FIG. 3. The (normalized) decay rate T(t), see Eq. ([TT}, of 
a massless scalar field (j> as it flows through the black hole's 
horizon. The units are given in terms of {GM)~^, which for 
a typical mass of W^Mq gives F ~ 0.2s~^. The runs were 
performed for different values of the amplitude A and of the 
mean wavenumber fco of the wave packet ((7|, whereas the 
width was fixed to ct = 5. The total interaction time between 
the black hole and the wave packet is around tr — 4a" in the 
two cases. The decay rate increases for larger values of fco, 
but it is insensitive to the values of the amplitude A, i.e., to 
values of the scalar field's strength. 



because of the possibility that SMBH could have grown 
because of the accretion a Quintessence-type scalar field. 
This is not the only possible case, but we can ask the 
same question about any other cosmological scalar field 
living around a black hole. 

We have explored the simplest possibility, that of a 
massless scalar field, for the motion of a wave packet in 
a fixed black hole spacetime using an EF gauge. To have 
a better visualization of the absorption of the scalar field 
by the black hole, we took advantage of the fact that 
one can write out a continuity equation from the KG 
equation. 

The corresponding conserved charge is the total mass 
of the wave packet, but more important is the defini- 
tion of a current density, with the help of which we were 



FIG. 4. Total absorption ratio A, see Eq. (fT^ . for different 
widths of the wave packets; for all runs fco = 0. The amplitude 
A was adjust in each case as to have always the same total 
mass M^ at the initial time. This is not strictly necessary, as 
the decay rate is independent of the packet's amplitude. We 
also show the fitting function p3p . 



able to define a (normalized) decay rate for the wave 
packet, whose magnitude is given by the black hole's 
mass, r(i) ~ {GM)~^. This means that less massive 
black holes accrete scalar field matter at a larger rate; 
for example, a black hole as massive as the Sun would 
accrete at an incredible rate of F ~ 10^s~^! In terms of 
the decay rate, we too found that the absorption depends 
on the mean wavenumber of the wave packet; actually, 
full absorption is reached for mean wavelengths smaller 
than the Schwarzschild radius, ko < rg . 

However, a new result showed up: the decay rate does 
not depend on the scalar field's strength. Moreover, we 
could use this result to show the dependence of the ab- 
sorption on the packet's width. By a limiting procedure 
on a fitting function, we determined the maximum to- 
tal absorption of a wave packet with a width much larger 
than the black hole's horizon: around e~°-^^ ~ 0.51. This 
is the result that may have relevance for cosmology, as 
we expect cosmological scalar fields to have very large 
intrinsic length scales (fco ~^ and cr — > oo) as compared 
to the Schwarzschild radius of supermassive black holes. 

In the massless case studied here, there were two length 
scales involved: the mean wavelength Xq = kQ , and the 
width of the wave packet a. We were able to show the 
general dependence of the (normalized) absorption rate 
on these length scales. In general, we can say that black 
holes are quite efficient in absorbing scalar fields, even in 
the case of very wide packets. 

The method outlined here can be extended to the mas- 
sive case. However, in the case the scalar field has a mass 
m, additional length scales appears in the problem in the 
form of the Compton length of the field, Xc = m~^ and 
the Schwarzschild radius rs (in the massless case, the 
Schwarzschild radius does not appear explicitly in the 



equations of motion) that may introduce non-trivial fea- 
tures in the motion of the wave packet and its absorption 
rate. This is ongoing research that we expect to report 
elsewhere. 
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